©Two Mathematical Models of Nonlinear Vibrations 

Model parameters are fit to empirical vibration data. 
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Two innovative mathematical models 
of nonlinear vibrations, and methods of 
applying them, have been conceived as 
byproducts of an effort to develop a 
Kalman filter for highly precise estima- 
tion of bending motions of a large truss 
structure deployed in outer space from a 
space-shuttle payload bay. These models 
are also applicable to modeling and 
analysis of vibrations in other engineer- 
ing disciplines, on Earth as well as in 
outer space. 

The first model is denoted the ampli- 
tude-dependent stiffness (ADS) model to 
emphasize the difference between it and 
the classical linear harmonic-oscillator 
model, in which stiffness is a constant. The 
ADS model is embodied in the equation 

x + £x + K{x,x)x = 0, 

where * is the instantaneous amplitude 
of the oscillating position or modal co- 
ordinate, ^ is a damping parameter, 
and K{x,x ) is the ADS. 

In the initial outer-space application, 
the ADS was represented by the follow- 
ing nonlinear function: 

K{x,x) = a + bA(x,x) + cA(x,x) , 

where a, b, and c are constant parame- 
ters to be obtained by fitting the model 
to empirical amplitude-versus-fre- 
quency data, and A(x,x) is a modal am- 
plitude. The amplitude-versus-fre- 
quency data are obtained by means of a 
moving-window estimation technique 
in which one analyzes the instanta- 
neous vibration waveform during a 
time window of about 90 percent of the 


time-average vibration period. The am- 
plitude and frequency are taken to be 
those of a sinusoid that makes the least- 
squares best fit to the instantaneous am- 
plitude during the window (see figure). 
The window is then moved by about 2 
percent of the average period and an- 
other best-fit sinusoid is found. This 
process is repeated until a suitably repre- 
sentative sample of the vibration wave- 
form has been acquired. 

The modal amplitude is given by 


A(x,x) 


x 2 + 


K(x,x) 


where K(x,x) is any reasonable approx- 



A Nonlinear Decaying Waveform is approximated 
with a best-fit sinusoid during a moving window. 
The resulting sinusoidal amplitude and frequency 
data are collected from fits for the entire sequence 
of window positions and used to characterize the 
frequency versus amplitude of the nonlinear wave- 
form. These frequency versus amplitude data are 
then fit to an amplitude-dependent stiffness (ADS) 
representation. 


imation of K{x,x) . One can refine the 
approximation iteratively, starting from 
K(x,x) = a, then using the resulting 
value of A(x,x) in computing a value of 
A'(x,x) by use of the above equation for 
A'(x,x ) . 

The second model, denoted the mo- 
ment-expansion (ME) model, is embod- 
ied in the equation 

x + M(x,x) = 0, 

where the function M(x,x) is a moment 
expansion that captures damping and 
stiffness effects. The moment expansion 
is given by 

3 3 

M(x,x) = y y s P J p ij x t P, 

j= 0 *=o 

where both i and j range from 0 to 3, 
except that there is no (i,j) = (0,0) 
term. In the original outer-space appli- 
cation, the parameters py are obtained 
from (1) modal position and velocity 
estimates obtained from Kalman-filter 
states and (2) derived accelerations. 

In a test relevant to the original 
outer-space application, the ADS and 
ME models were compared with each 
other, with a linear model, and with a 
prior nonlinear model known as the 
Duffing model. The ADS model was 
found to yield the least error. 

This work was done by Paul Brugarolas, 
David Bayard, John Spanos, and William 
Breckenridge of Caltech for NASA’s Jet 
Propulsion Laboratory. For further informa- 
tion, contact iaoffice@jpl.nasa.gov. NPO- 
41360 


@ Simpler Adaptive Selection of Golomb Power-of-Two Codes 

The selected code-parameter value is within 1 of the 
optimum value. 
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An alternative method of adaptive se- 
lection of Golomb power-of-two (GP02) 
codes has been devised for use in effi- 
cient, lossless encoding of sequences of 
non-negative integers from discrete 
sources. The method is intended espe- 
cially for use in compression of digital 
image data. This method is somewhat 
suboptimal, but offers the advantage in 


that it involves significantly less compu- 
tation than does a prior method of adap- 
tive selection of optimum codes through 
“brute force” application of all code op- 
tions to every block of samples. 

A rather lengthy discussion of back- 
ground is necessary to give meaning to a 
brief summary of this innovation. For 
positive integer, m, the rath Golomb 


code defines a reversible, prefix-free 
mapping of non-negative integers to 
variable-length binary code words. 
Golomb codes are optimum for geomet- 
rically distributed sources (a model that 
frequently arises in image compression) : 
In the case of a geometrically distributed 
random variable, 8, the appropriately se- 
lected Golomb code minimizes the ex- 
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pected code-word length over all possi- 
ble lossless binary codes for 8. 

In a GP02 code, m = 2 k , where A is a 
non-negative integer. Such a code makes 
the coding process particularly simple: 
The code word for the integer 8 consists 
of the unary representation of 1. 8/2*1 (that 
is, 18/2*1 zeros followed by a one) concate- 
nated with the k least significant bits of 
the binary representation of 8. More 
specifically, the code is called a GP02 
code of parameter k. 

The problem is to calculate or esti- 
mate the value of code parameter k that 
minimizes the expected bit rate (the av- 
erage number of encoded bits per 
source symbol) for an image or other 
source. This problem arises in Rice cod- 
ing, which is a coding method well 
known among experts in data compres- 
sion. The Rice algorithm encodes a 
block of samples by use of the best code 
option for the block from among several 
candidate codes that consist mosdy of 
different GP02 codes. A fixed number 
of bits are used preceding the encoded 
block to indicate which code was se- 
lected. The Rice method does not spec- 



Upper and Lower Bounds on the optimum value 
of the code parameter make it possible to re- 
duce the number of code options that must be 
considered. 

ify how to find the best code option, and 
the most common approach is to ex- 
haustively try every code option to pick 
the best one for each block. Information 
from previously coded blocks is not uti- 
lized. This concludes the background in- 
formation. 

In the present method, unlike in the 
Rice method, one utilizes the mean 
sample value in each block. The 
method is based partly on a theoretical 


derivation of bounds on the optimum 
value of k as functions of the mean 
sample value (see figure). These 
bounds are such that no more than 
three code choices can be optimum 
for a given mean sample value. For a 
given mean value, one of the three 
candidate codes is selected in a proce- 
dure that involves only integer arith- 
metic (without divisions) and table 
look-ups. It has been shown that the 
value of k selected in this relatively sim- 
ple procedure is always within 1 of the 
optimum k value for the source, and 
that the cost added by the suboptimal- 
ity of the selection is never more than 
1/2 bit per sample and no more than 
about 13-percent inefficiency. In prac- 
tical image compression experiments, 
the cost added by the suboptimality of 
the selection is negligible. 

This work was done by Aaron Kiely of Cal- 
tech for NASA’s Jet Propulsion Laboratory. 

The softzvare used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-41336. 
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